Non-Markovianity of local dephasing channels and time invariant discord 
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We study non-Markovianity and information flow for qubits experiencing local dephasing with an 
Ohmic class spectrum. We demonstrate the existence of a temperature-dependent critical value of 
the Ohmicity parameter s for the onset of non-Markovianity and give a physical interpretation of 
this phenomenon by linking it to the form of the reservoir spectrum. We demonstrate that this link 
holds also for more general spectra. We unveil a class of initial states for which discord is forever 
frozen at a positive value. We connect time invariant discord to non-Markovianity and propose a 
physical system in which it could be observed. 
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Introduction.- Qubits subjected to local purely dephas- 
ing noise are ubiquitous models of open quantum sys- 
tems, and they have been studied extensively in the liter- 
ature. Examples include dephasing in quantum registers 
1t3(j ultracold gases [J,|5|, quantum metrology protocols 
6|, quantum biological systems |7|, and dynamical de- 
coupling theory [8[ . From a theoreticalpoint of view the 
dephasing model is exactly solvable [l|-|3[ , and hence it is 
an ideal testbed to investigate one of the most thrilling 
fields of the theory of open quantum systems, that of 
non-Markovian quantum processes |9| . 

Recently, a great deal of attention has been devoted 
to the study of systems whose reduced dynamics are 
characterised by memory effects and recoherence phe- 
nomena, emerging from non-trivial correlations with 
an environment. Such dynamics are typically called 
non-Markovian. Memory effects and non-Markovianity 
have been shown to be a resource for quantum tech- 
nologies [y, Ild - ll3{ and consequently measures of non- 
Markovianity have become important as quantifiers of 
this resource [l4Ml6[ | . Moreover, it has been shown that 
non-Markovianity of a quantum probe can indicate a 
quantum phase transition occurring in a complex envi- 
ronment, with which the probe is interacting [17] . 

Non-Markovian features play an important role in sys- 
tems where the frequency spectrum of the environment 
is structured. However, a connection between the gen- 
eral form of the spectrum and the memory effects in the 
reduced system dynamics has not been elucidated un- 
til now. In this Letter we establish this connection by 
unveiling a necessary condition on the form of the spec- 
trum to induce non-Markovian dynamics for a dephasing 
qubit. We then focus on the widely used Ohmic class of 
reservoir spectra and show that the condition is both nec- 
essary and sufficient for this type of spectra. Moreover, 
we demonstrate that only super-Ohmic environments can 
induce non-Markovian dynamics. This means that even 
if the reduced dynamics is exact, and hence no Markovian 
approximation has been performed, the time evolution of 
the qubit does not present any memory effects or reco- 



herence for Ohmic and sub-Ohmic spectra. 

Non-Markoivan dynamics can prolong the existence of 
quantum properties, thus delaying the quantum to clas- 
sical transition. We present here a striking example of 
this phenomenon by studying the decay of quantum and 
classical correlations in a non-Markovian open quantum 
system. When two qubits interact with purely dephasing 
local environments, the dynamics of both classical and 
cmantum correlations can exhibit sudden changes [19l — 
21] ; for certain initial conditions one observes a sudden 
transition from classical decoherence (decay of classical 
correlations) to quantum decoherence (decay of quantum 
correlations), which is preceded by a finite interval of 
time when quantum discord, a commonly used measure 
of quantum correlations 22|, |23| , is frozen to a non-zero 
value [20( . This scenario holds for a Markovian model of 
noise, that is, when one describes the system by means 
of a master equation of Lindblad-Gorini-Kossakowski- 
Sudarshan form 18|, and it has also been observed for 
non-Markovian random telegraph noise [24j and several 
other physical models [251 ] - 

It is natural to wonder how this transition behaves 
in the exact pure dephasing model considered here, and 
ask what the role of non-Markovianity is in this process. 
In this Letter we demonstrate that for this model the 
sudden transition does not always occur and the quantum 
correlations can behave in a time invariant way, that is, 
remaining constant at all times, while the state of the 
system and all other dynamical quantities evolve. We 
present the conditions for permanently frozen discord and 
discuss its microscopic origin. Specifically we point out 
how time invariant discord is related to non-Markovian 
features of the exact dynamical map. 

Non-Markovianity.- Let us consider the following mi- 
croscopic Hamiltonians describing the local interaction of 
a qubit and a bosonic reservoir, in units of h, 

H = lo o- z + ]P uj k a\a k + ^ o~z{gkak + 9k a t)> 

k k 

with wq the qubit frequency, ujk the frequencies of the 



reservoir modes, a^ (a k ) the annihilation (creation) oper- 
ators and gk the coupling constant between each reservoir 
mode and the qubit. In the continuum limit J^k \9k\ 2 ~^ 
J dujJ(uj)S(u!k — w), where J(uj) is the reservoir spectral 
density. This model can be solved exactly [lrlsl- The 
master equation for the qubit, in the interaction picture, 
is given by 



p = j{t)[a z pa z - p]/2, 



(1) 



and is time-local. If the environment is initially in a 
thermal state, the time-dependent dephasing rate takes 
the form 

-f(t) = / dwJ(w) coth [hbj/2k B T] sin(wt)/u;, (2) 
resulting in the decay of the density matrix off-diagonal 



elements: Pij{t) 



-A(t] 



Pij(0), i ^ j, with dephasing 



factor A(t) = 2 £ j(t')dt' given by 
A(t) = 2 / du> J(w) coth [huj/2k B T] [1 - cos(ujt)} / uj 2 



dw g(w, T) [1 — cos(wi)] 



(3) 



For this model the non-Markovianity measures based 
on information flow F|, on entanglement with an an- 



cilla [15| and on Fisher information [16[ all predict the 



same crossover between Markovian and non-Markovian 
dynamics. The crossover is signaled by the onset of peri- 
ods during which the dephasing rate is negative or, equiv- 
alently, the dephasing factor A(i) always decreases with 
time. 

To elucidate the origin and physical meaning of non- 
Markovianity in dephasing channels we notice that, fol- 
lowing Ref . [2| , one can describe the effect of the qubit 
on its environment by a displacement operator acting 
on each environment mode, with the associated phase 
conditional on the state of the qubit. The two qubit 
states excite each mode with opposing phases, leading to 
a decrease in the overlap between the states of the mode 
in each case; this is the physical cause of decoherence. 
Destructive interference between excitations of a mode 
at different times then leads to rccohercnces at the fre- 
quency of the mode; it is the balance between these two 
effects for different modes, captured exactly by Eq. ([3]), 
that determines whether the dynamics is non-Markovian. 

Using Eq. (J3j> , we establish a simple link between the 
onset of non-Markovianity and the form of the reservoir 
spectrum. As the cosine transform of a convex function 
is monotonically increasing, we deduce that a sufficient 
condition for Markovianity is that g(u,T) is convex or, 
equivalently, the non-convexity of g{ui, T) is a necessary 
condition for non-Markovianity. Physically, a convex 
g{uj,T) means that any recoherence is always outweighed 
by more decoherence from lower frequency modes. As we 



now show, this condition helps classify the dynamics for 
Ohmic-likc spectral densities of the form 



J(uj) 



o/u 



(4) 



where uj c is the reservoir cutoff frequency. By chang- 
ing the s-parameter one goes from sub-Ohmic reservoirs 
(s < I) to Ohmic (s = I) and super-Ohmic (s > 1) reser- 
voirs, respectively. We stress that such engineering of 
the Ohmicity of the spectrum is possible when simulat- 
ing the dephasing model in trapped ultracold atoms, as 
demonstrated in Ref. [5[. A closed analytic expression 
for the time-dependent dephasing rate can be found in 
both the zero T and the high T limit. In the former case 
one obtains 

7o(t, s) = [1 + (uj c t) 2 }- s/2 T[s} sin [s arctan(o; c t)] , (5) 

with r[x] the Euler gamma function. For high T, instead, 
jHT(t, s) = 2fc B T7o(t, s - 1)/lo c . 

Starting from Eq. ([5]) it is straightforward to prove 
that at zero T the dephasing rate takes temporarily neg- 
ative values if and only if s > s cr ,;t = 2. Hence, mem- 
ory effects leading to information back flow and recoher- 
ence occur only if the reservoir spectrum is super-Ohmic 
with s > 2. Equally, Eq. (0 leads to s cr i t = 3 for 
high T. Moreover, we have established numerically that 
s cr it increases monotonically with the temperature until 
it reaches its maximum value s cr i t = 3 at infinite temper- 
ature. The existence of a temperature dependent critical 
value of the Ohmicity parameter, ruling the Markovian 
to non-Markovian transition, is one of the main results 
of this Letter. We now explain this result in terms of the 
reservoir spectrum. 

It can be shown analytically that the integrand of Eq. 
Q for the Ohmic class becomes a non-convex function 
of uj for s > s cr it for both zero and infinite T. For inter- 
mediate temperatures, numerical investigation also indi- 
cates that the value of s for which the function g(u>,T) 
changes from convex to non-convex coincides with s cr u, 
implying that the condition on the non-convexity of the 
spectrum is necessary and sufficient for non-Markovianity 
at all T. Therefore, not only does convexity guarantee 
decoherence always outweighs recoherence, but for these 
systems, it is required; this highlights the key role of the 
low frequency part of the spectrum in the occurrence of 
information backflow. 

Time invariant discord.- For a bipartite state pab the 
quantum discord Q(pab) is defined as the difference be- 
tween the total correlations of the system, given by the 
mutual information I{pab) = S(pa) + S(pb) — S(pab) 
with S(p) = — Tr(plogp) the von Neumann entropy, and 
the classical correlations C(pab) = max.n A [J(HApAB)} 
|22l [23j. In the latter expression the maximization is 
performed over all sets of orthogonal projections 11^ on 
state A and J describes the effect of the projective mea- 



surcments on the state. In general the maximization pro- 
cedure makes calculations of the quantum and classical 
correlations difficult, however in this work we focus on a 
class of states for which the maximizing measurement is 
known and quantum discord has an analytical expression 

In Ref. [20j it was shown that if two qubits interact 
locally with Markovian dephasing environments, there 
exist classes of states for which during an initial time in- 
terval < t < t quantum discord remains constant while 
classical correlations decay. For t > t, on the other hand, 
quantum discord decays while classical correlations re- 
main constant. In some sense this is counter- intuitive as 
we would expect quantum properties to start decaying 
before classical properties, since they are the most sen- 
sitive to the deleterious effects of the environment. This 
phenomenon, named the sudden transition from classical 
to quantum decoherence, was investigated originally the- 
oretically and experimentally when the qubit dynamics 
is described by a Markovian master equation of the form 

PA(B) = l[<rf PA{B)<rf ~ PA(B)]f2, wit h 7 > the 
constant dephasing coefficients, a- the Pauli operator 
in direction j acting on qubit A(B), and j = x, y, z. For 
the sake of concreteness we consider dephasing along the 
z- direction only. The sudden transition occurs for initial 
Bell diagonal states of the form 



Pab = 



(1 



|*±)(tf= 



(1-c) 



4>±)($ d 



(6) 



where |#±) = (|00) ± |11))/V2 and |$±) = (|01) ± 
|10))/-\/2 are the four Bell states and \c\ < 1. The sudden 
transition time is given by i = — ln(|c|)/(27). Therefore, 
it is possible to increase the time interval over which the 
discord is constant. For increasing values of t, however, 
the discord decreases towards its zero value obtained for 
c = 0. 

We now study whether the sudden transition from clas- 
sical to quantum decoherence occurs also for the exact 
model of dephasing considered here and, if so, the effect 
of the reservoir spectrum on the behavior of quantum and 
classical correlations. For the initial class of states of Eq. 
(J6j) the mutual information and classical correlations take 
the form 



Z[pAB(t)} 



l + (-lV'c 



2 

E 

3=1 

y i + (-i)V 

^ 2 

3=1 



log 2 [l + (-l)'c] 



(7) 



•A(t) 



and 



C[pAB(t)} 



2 

E 

3=1 



i + (-i) j x(t) 



log 2 [l + (-l)^e- A W] 



log a [l + (-l)'x(*)],(8) 



where %(£) = max{e A '*\c}, and we have taken c posi- 
tive for the sake of simplicity. 
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FIG. 1. (Color online) (a) Landscape of correlation dynamics 
in the s — f-plane, for c = 0.1 and T = 0. Blue areas denote 
parameters (£, s) corresponding to classical decoherence, red 
areas to quantum decoherence and the intersection between 
the two, marking the values of s and t satisfying Eq. Q, 
defines the transition time (as a function of the reservoir 
spectrum parameter s. The two insets show discord (solid 
red line) and classical correlations (dashed black line) for two 
specific choices of s; (b) for s — 1 the system has a sudden 
transition from quantum to classical decoherence while for (c) 
s = 2.5 the discord is frozen forever. The blue dot in (a) and 
(b) points the transition time i for s = 1. 



From Eqs. (O-© one sees immediately that when 
e -A({) ^ c |- nc c i ass i ca [ correlations decay while the dis- 
cord, given by the first term in Eq. (|7|), remains constant. 
On the other hand, if a finite transition time t such that 



-Hi) 



(9) 



exists, then for t > i the discord starts decaying and 
the classical correlations stay constant. Contrary to the 
Markovian dephasing model the transition time t now 
crucially depends not only on c but also on the parameter 
s and on the reservoir temperature T through A(i). 

Figure 1 shows the values of s and t for which condi- 
tion ^ is satisfied, for c = 0.1 and T = 0. For a certain 
range of the parameter s Eq. © has a solution and ac- 
cordingly the system has a sudden transition from classi- 
cal to quantum decoherence at time t. Interestingly, we 
also discover a range of values of s for which Eq. © has 
no solution and the transition time t does not exist. For 
these values of s only classical correlations are affected 
by noise, leading to classical decoherence, while discord 
remains frozen forever. The two different cases are illus- 
trated in Figs. 1(b) and (c) where we plot the classical 
correlations and discord for the Ohmic case s = 1 and 
s = 2.5, respectively. 

For the zero T case, the possibility of having time in- 
variant discord depends on both the initial state of the 
two-qubit system, that is on the parameter c, and on 
the parameter s, describing the structure of the reservoir 
spectral density. By looking at the asymptotic long time 
limit of condition ^ we can define the s and c param- 
eter space for which time invariant discord exists. This 
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FIG. 2. (Color online) The shaded region marks range of 
parameters s and c for which the discord is frozen forever 
for T — 0. Outside this region one will always observe a 
transition from classical to quantum decoherence. 



is shown in Fig. 2. Note that the value of the frozen 
discord is Q= (1 + c) log 2 (l + c)/2 + (l-c) log 2 (l-c)/2, 
which is very small for small values of c, hence it can be 
argued that the phenomenon of time invariant discord is 
significant for larger values of c, roughly corresponding 
to 2 < s < 3, for example, when the dynamics is non- 
Markovian. Increasing the temperature rapidly destroys 
the time invariant discord phenomenon. In the high tem- 
perature limit, indeed, one cannot find any value of s and 
c for which this effect occurs. We show in the following, 
however, that there exist realistic physical settings for 
which the effect can be observed. 

Non-Markovianity and time invariant discord. - 
Let us now link the occurrence of time invariant dis- 
cord in the two qubits dynamics to the form of the reser- 
voir spectrum and non-Markovianity in the single qubit 
dynamics. It is easy to convince oneself that the time 
invariant discord phenomenon can occur only for reser- 
voir spectra leading to a bounded value of A(i). This 
ensures the existence of values of c such that e w > c 
for all t, implying that Eq. © is never satisfied. An 
asymptotic divergence of A(i), on the contrary, allows 
for the existence of a transition time t. Such a diver- 
gence, and therefore absence of time invariant discord, 
rests on the divergence of g(u>, T)/u when uj — > occur- 
ring for s < 1(2) at zero (finite) temperature. Similarly, 
convexity and thus Markovianity is ensured if g(u, T) di- 
verges at low frequencies, occurring for s < s crit = 2(3). 
The above demonstrates that time invariant discord and 
non-Markovianity are intimately related and ultimately 
rely on the eventual dominance of recoherence over de- 
coherence; thus both require the suppression of coupling 
to low frequency modes, embodied by the low frequency 
dependences of J(w) and g(w,T). 

Finally, we note that a possible experimental setup in 
which the time invariant discord could be observed is 
an array of double-well impurities immersed in a Bosc- 
Einstcin condensate, as discussed in Rcfs. [4|, |5j. When 



the impurities are far apart, that is, the distance D be- 
tween the impurities is much greater than the distance 
L between the two potential wells forming an impurity, 
the collisions with the ultracold gas lead to an effec- 
tive local pure dephasing model as the one considered 
in this Letter. For the parameters of Fig. 5 in Ref. [4| , 
with D = 20L, the maximum value reached by the de- 
coherence factor before attaining its stationary value is 
max t A(£) ~ 0.058, hence time invariant discord occurs 
for < c < 0.94, that is for a wide range of initial states. 
As the discord increases for increasing values of c, within 
this physical systems one can freeze discord to values 
close to its maximum forever. Moreover, these systems 
have proven to be very resistant to the effects of finite 
temperature [23], and we therefore expect to find high 
values of discord also for realistic temperatures of the 
order of T = 10 - lOOnK. 

Conclusions.- Non-Markovianity reflects the ability of 
an open system to regain and retain quantumness, since 
previously lost quantum information can partly flow back 
into the system and, in some cases, be trapped. As reser- 
voir engineering techniques become experimentally feasi- 
ble, it is crucial to establish qualitative and quantitative 
links between the occurrence of non-Markovianity and 
the form of the environmental spectrum. No such con- 
nection was known until now. In this Letter we have 
paved the way to these studies by presenting a neces- 
sary condition on the form of the spectrum for the non- 
Markovianity of a qubit undergoing pure dephasing. We 
have proven that, for the Ohmic-class of spectral den- 
sities the condition is necessary and sufficient, and dis- 
covered the existence of a temperature-dependent criti- 
cal value of the Ohmicity parameter for the Markovian 
to non-Markovian crossover. For two qubits in locally 
dephasing environments, we have unveiled a new phys- 
ical phenomenon, time invariant discord, and explored 
its relation to the non-Markovianity of the individual de- 
phasing qubits. 

Our results shed light on the physical origin of de- 
coherence in a paradigmatic open quantum systems 
model, and are therefore of great fundamental impor- 
tance. Moreover, as the model describes the dominant 
source of noise in several systems used for quantum tech- 
nologies, our results may lead to the implementation of 
more resistant quantum protocols and devices. 
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